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By  a slight  modification  of  Saffman-Taylor’s  viscous  finger,  we  remove  its  singu- 
larity and  present  a finite-length  nonsingular  Saffman-Taylor  finger.  Its  dynamic 
evolution  solutions,  both  one-step  theory  prediction  and  numerical  simulation  in 
the  long-time  limit  with  zero  surface  tension,  give  an  alternative  proof  for  the  an- 
alytical selection  demonstration  of  the  Saffman-Taylor  finger  width  in  the  absence 
of  surface  tension.  As  a simple  example  of  pattern  selection,  which  has  real  exper- 
imental background,  this  work  not  only  contradicts  the  generally  accepted  belief 
that  surface  tension  is  indispensable  for  the  selection  of  the  t-width  finger  but  also 
provides  more  models  to  compute  the  evolution,  competition  and  ramification  of 
multiple  fingers  numerically  in  straight  channel  as  well  as  circular  disc  geometry. 


1 Introduction 

The  Saffman-Taylor  problem  1 has  played  a central  role  in  the  study  of  viscous 
fingering  in  a Hele-Shaw  cell  2,  which  was  modeled  by  two-dimensional  potential 
flow  at  the  interface  between  two  fluids.  Saffman  and  Taylor  found  analytically  a 
continuous  family  of  steady-state  solutions,  which  shows  fingers  of  different  width 
could  exist  in  the  absence  of  surface  tension,  but  their  experiments  with  negligible 
surface-tension  effect  in  1958  and  numerical  calculations  made  by  McLean,  Saffman 
and  Vanden-Broeck  in  the  presence  of  surface-tension  effect  in  1981  and  1983  showed 
no  hint  of  the  continuous  family  3,4 . On  the  contrary,  the  finger  with  relative  width 
A = 0.5  is  always  singled  out  at  the  zero  surface  tension  limit.  This  problem 
(which  is  so  called  Saffman-Taylor)  has  been  of  much  subsequent  interest  because 
it  is  universal,  i.e.,  the  same  selection  phenomenon  is  common  for  displacement  of 
various  viscous  liquids  by  less  viscous  ones  for  immiscible  incompressible  liquids. 
Much  work  5 was  done  toward  solving  this  Saffman-Taylor  puzzle.  However,  it 
has  been  widely  accepted  that  the  inclusion  of  surface  tension  is  the  only  way 
to  select  the  most  stable  finger  width  since  Saffman-Taylor  proposed  that  surface 
tension  between  the  two  fluids  would  solve  the  selection  problem  1 . Although  several 
works  6 confirmed  numerical  evidence  of  the  discrete  spectrum  of  A,  decreasing  to  A, 
in  the  limit  of  low  surface  tension,  exact  solutions  are  rare  with  zero  surface  tension. 
Recently,  Mineev- Weinstein  7 analytically  solved  the  finger  selection  problem  in  the 
absence  of  surface  tension.  He  showed  an  exact  result  that  a generic  interface  in 
a Hele-shaw  cell  evolves  to  nonlinearly  stable  single  uniformly  advancing  finger 
occupying  one-half  of  the  channel  width,  which  contradicts  the  generally  accepted 
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Figure  1. 


belief  that  surface  tension  is  indispensable  for  the  selection  of  the  |-width  finger. 

We  find  that  the  problem  of  the  finger  width  selection  stems  from  the  singularity 
of  Saffman- Taylor  finger,  which  is  half  infinite  in  length.  The  singularity  makes  the 
numerical  calculation  impossible  and  is  not  in  agreement  with  the  experimental 
results.  We  introduce  a positive  variable  e(t)  to  remove  the  singularity.  By  the 
conformal  mapping  method  presented  by  Bensimon  8 , we  predict  that  there  exists 
a unique  steady  width  A = 0.5  of  the  viscous  finger  to  which  the  fingers  of  other 
width  approach  in  the  absence  of  surface  tension  after  a one-step  analysis.  Then  we 
verify  this  prediction  by  numerical  computation  and  show  its  long  time  evolution. 


2 Nonsingular  Saffman- Taylor  Finger  and  One-step  Dynamic 

Evolution  Theory 

Following  Saffman  and  Taylor  we  consider  the  analysis  of  a long  air  bubble  or 
finger  moving  through  a channel  in  the  Hele-Shaw  cell  bounded  by  two  straight 
parallel  walls.  We  take  the  y-axis  (Fig.  1)  perpendicular  to  the  walls  and  the  x-axis 
along  the  central  axis  of  symmetry  of  the  channel.  The  separation  of  the  walls  is 
normalized  to  27r  and  the  velocity  of  the  fluid  at  infinity  in  front  of  the  finger  is 
normalized  to  unity.  Let  z = x + iy  and  $ = <j>  + iip,  <f>(z)  and  ip(z)  being  the 
velocity  potential  and  stream  function,  respectively. 

In  case  of  the  absence  of  surface  tension,  the  interface  is  described  by  0 = 0, 
which  can  be  conformally  mapped  to  the  circumference  of  a unit  disk  in  the  complex 
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plane  ( = e *.  By  these  notations  the  Saffman- Taylor  solution  1 can  be  written  as 

z = 2(1  - A)ln^(e~t1p  + 1)  + iip,  (1) 

where  the  single  parameter  A depicts  the  width  of  the  finger.  Saffman-Taylor  finger 
is  half  infinite  in  length  and  the  solution  is  steady  in  the  frame  co-moving  with  the 
finger  tip  which  has  a speed  1 relative  to  the  laboratory  frame. 

In  order  to  remove  the  singularity  at  ip  = ±7r  in  Eq.  (1),  we  choose  the  solution 
in  the  following  form  in  the  laboratory  frame: 

z = 2(1  - A)2«i(C  + 1 + e)  - lnC+  j + z((,t),  (2) 

where  e = e(t)(=e o at  f=0)  is  a small  positive  real  number  varying  with  time  f, 
z(C,t)(— 0 at  t=0)  describes  the  possible  additional  time-dependence  of  the  finger. 
The  functional  form  of  e(t)  and  z(£,  t)  are  to  be  determined  by  the  evolution  equa- 
tion of  interface.  In  Eq.  (2)  the  singularity  of  Saffman-Taylor  solution  at  £ = — 1 
on  the  unit  circle  is  moved  to  outside  of  the  unit  disk  and  therefore  the  length  of 
the  finger  becomes  finite  (Fig.l). 

The  corresponding  equation  of  the  boundary  of  the  nonsingular  Saffman-Taylor 
finger  is 


z = 2(1  - A )ln^(e  ^ + 1 + e)  + iip  + ~ + z(ip,  t)  (3) 

At  moment  f , the  finger  tip  is  located  at  ip  = 0,  (j>  = 0 and  the  rear  end  at  ip  — 
±7r,  (p  = 0,  therefore  the  length  of  the  finger  is 

L(t)  = £0  — aq  = 2(1  — X)ln(l  + -)  (4) 


In  our  notation  Bensimon’s  equation  8 in  case  of  the  absence  of  surface  tension  is 
written  as 


dz 

dt 


(5) 


where 


G(C)  = A{g(ip)},  g{ip)  = Or-  (6) 

and  the  suffix  T means  that  the  value  of  the  function  in  the  square  parentheses  is  to 
be  taken  on  the  unit  circle  where  Q = , A is  an  operator  analytically  continuing 

the  real  function  g (ip)  defined  on  the  circumference  of  the  unit  circle  as  the  real 
part  of  complex  function  G(£)  in  the  interior  of  it.  The  analytical  continuation  is 
performed  by  Fourier  expansion. 

When  t= 0,  Eq.  (2)  becomes 

z = 2(1  - A)Zn-(C  + 1 + eo)  - (7.1) 

and  the  equation  of  the  nonsingular  Saffman-Taylor  finger  is 
z = 2(1  - A )Zn^(e-1^  + 1 + e0)  + iip 


(7.2) 


Then 


, _ 1 + 2(1  + £o )cosip  + (1  + eo)2 

9[V)  ~ (1  - 2A)2  - 2(1  - 2A)(1  + eo)cos1>  + (1  + c0)2 
taking  a trial  solution  as 

rtr)  = a + bl' 

(l-2A)C-(l  + eo) 

with  a,  6 to  be  decided  , we  have 

1?  rfrll  ^ + a)(!  ~ 2A)  “ b(l  + £o )cosip  - a(  1 + e0) 
ig|r  ~ (1  - 2A)2  - 2(1  - 2A)(1  + e0 W + (1  + e0)2 ' 

Equating 

ReG(i/))\r  =g{4>) 


we  get 

(1  — 2A)  + (1  — 2A)(1  + eo)2  + 2(1  + eo)2 
“ (1  — 2A)2  — (1  + e0)2 


(8) 

(9) 


(10) 

(11) 

(12) 


_ (2(1  — 2A)  + (1  + eo)2  + 1)(1  + £o) 

(1  - 2A)2  - (1  + e0)2 

Substituting  eqs.(12)  and  (13)  in  G(( '),  we  get 

_r— rtrl  I _ (l  - A)(i  + c0)[i  - (l  + e0)2] 

A A[(l-2A)2-(l+eo)2](C+l  + eo)' 

(l-2A)(l-A)[l-(l  + £o)2]C 
A[(l-2A)2-(l  + e0)2](C+l  + e0)' 
Meanwhile,  from  Eq.(2)  at  time  t,  we  have 

dz  1 2(1  — A)  deo  dz 

dt  A~^C  + 1 + e dt  ^ dt 

Comparing  Eq.(14)  and  Eq.(15),  we  have 

2(1  — A)  de o _ (1  — A)(l  + eo)[l  ~ (1  + eo)2] 

C + 1 T Co  dt  A[(l  — 2A)2  — (1  + £o)2](C  T 1 + cq) 


(13) 


(14) 


(15) 


(16) 


dz  _ (1  — 2A)(1  — A) [1  — (1  + c0)2]C 

dt  A[(l  — 2A)2  — (1  + eo)2](C  + 1 + eo) 
and  at  the  boundary 

dz_,  _ (1  - 2A)(1  - A)[l  - (1  + c0)2] { [1  + (1  + £q )cosip]  - i(  1 + e0)simp} 

dt  r A[(l  - 2A)2  - (1  + e0)2][l  + 2(1  + e0)cosip  + (1  + e0)2] 


In  the  first  time  step  we  have 


(18) 
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thus 

A9  = TiP~ At  = (1  - 2A)(1  - A)[l  - (1  + 6q)2][1  + (1  + £p )cosjj} 

dt  A[(l  — 2A)2  — (1  + e0)2][l  + 2(1  + eo)cosip  + (1  + e0)2] 


(20) 


= dz  (1  - 2A)(1  - A)[l  - (1  + ep)2](l  + e0)simp 

V dt  A[(l  - 2A)2  - (1  + e0)2][l  + 2(1  + eo)cosip  + (1  + e0)2] 


Choosing  a suitable  £q  from  Eqs.(20)  and  (21),  we  can  see 


(21) 


• If  originally  A > 0.5,  then  Ax  > 0,  Ay  < 0,  the  finger  is  being  stretched  to  be 
longer  and  thinner  in  later  time,  i.e.,  A tends  to  decrease. 

• If  originally  A < 0.5,  then  Ax  < 0,  Ay  > 0,  the  finger  is  being  compressed  to 
be  shorter  and  wider  in  later  time,  i.e.,  A tends  to  increase. 

• If  originally  A = 0.5, then  Ax  = 0,  Ay  = 0,  the  finger  keeps  its  width  forever. 


Thus  in  our  one-step  dynamic  evolution  theory  about  the  evolution  of  nonsingular 
Saffman-Taylor  finger,  we  conclude:  There  exists  a unique  steady  width  A = 0.5 
of  the  viscous  finger  to  which  the  fingers  of  other  width  approach.  This  is  exactly 
what  Saffman  and  Taylor  observed  in  the  experiment  in  1958. 


3 Numerical  Computation 


From  the  one-step  analysis  described  above  we  can  not  predict  completely  the 
shape  change  in  the  long  run,  however,  the  evolution  may  be  obtained  by  numerical 
computation.  At  time  t,  if  z(ip,  t ) = x(ij),  t)+iy(ip,  t ) is  given,  considering  Q = e~^ 
we  have 


[£<%z]r  = id^z 

substitution  of  eq.(22)  in  eq.(6)  yields 

9W  = (9*y)2  + (cV*)2 

According  to  the  Poisson  integral  formula—,  if 

CO 

y(V>)  = a0  + ^(a„ei^  + <e-i^) 

n—1 

we  have 

OO 

G(C)=A{y(r/0}=oo  + 2 5>nCn) 

n—1 


(22) 

(23) 


(24) 


(25) 


Noticing  the  symmetry  of  the  finger  and  substituting  Eqs.(22)  and  (23)  into  Ben- 
simon  equation  (5),  we  get  the  nonlinear  partial  differential  equation 

OO 

dtz  ~ -[2  y^(a nsinmp)  + ig(ip)]d^z 

n=l 


(26) 
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Figure  2.  Stable  evolution  of  NST  with  different  A's,e o = 0.3 


where 


„ _i  r 

Q>n  — / 

7T  Jo 


s(V') 


co snip  dip 


If  the  initial  condition  z(ip,  0)  = x(ip,  0)  + iy{ip,  0)  is  given,  we  can  compute  the 
equation  (26)  with  various  A and  get  a family  of  evolution.  Fig.  2 shows  the  evolu- 
tion of  nonsingular  Saffman-Taylor  finger.  It  can  be  seen  from  the  figure  that  the 
above  conclusions  are  valid  in  the  long  run  of  evolution. 

4 Conclusions  and  Prospect 

Using  the  conformal  mapping  method  presented  by  Bensimon,  we  present  a non- 

singular Saffman-Taylor  finger.  Both  theoretical  and  numerical  solutions  verify  the 

pattern  selection:  the  tracing  finger  always  tends  to  |-width  finger  in  the  long  run 
with  zero  absence  surface  tension  which  resolves  the  Saffman- Taylor’s  puzzle.  Being 
a simple  example  of  pattern  selection  which  has  real  experimental  background,  this 
work  contradicts  the  generally  accepted  belief  that  surface  tension  is  indispensable 
for  the  selection  of  the  width  finger,  a result  in  parallel  to  paper  in  ref 4 * *  7 . 
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Figure  3.  Straight  multiple  fingers  (left)  and  round  fingers  (right) 


Futhermore,  the  solution  (2)  of  single  finite  finger  can  be  generalized  to  cases  of 
multiple  fingers 

Z = [c-»*+i(n-l)ir  + l + e]  + $ (28) 

n 2 

and  of  round  fingers 

Inz  = + 1 + el  + $ (29) 

n 2 

(see  Figs.  3 and  4 respectively),  may  produce  more  interesting  pattern  selection 
problems  of  evolution,  competition  and  ramification  of  multiple  fingers  in  straight 
channel  as  well  as  circular  disc  geometry. 

It  should  be  mentioned  that  Bensimon’s  equation  is  a highly  nonlinear  equation 
which  meets  computational  difficulty.  Tiny  round-off  errors,  especially  the  residual 
errors  in  summing  up  the  Fourier  series  of  Eq.  (26)  will  be  magnified  quickly.  Special 
efforts  has  been  made  to  reduce  these  errors  to  extremely  low  level  in  order  to  make 
the  computation  possible,  which  has  been  completed  10 . If  surface  tension  exists,  we 
can  compute  Bensimon’s  equation  in  the  same  way,  but  special  method  of  numerical 
computation  must  be  designed. 

Pattern  selections  of  nonsingular  Saffman-Taylor  finger  in  the  presence  of  surface 
tension  and  multi-fingers  constitute  our  further  work. 
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